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Abstract 

The development of efficient protocols for pure and mixed states prepa¬ 
ration is challenging task. Most of the theory of quantum information 
applications has been developed for fully coherent or completely incoher¬ 
ent light. However, in many situations of interest partially coherent light 
has been proved to be a more robust model of radiation. In this paper 
the underpinning theory of two photon amplitude functions for down- 
converted fields with partially coherent pump beams is investigated. By 
using the generalised Siegert relations and the coherent mode representa¬ 
tion of the cross spectral density matrix the two photon amplitude is fully 
characterised for partially coherent beams. A number of correlation prop¬ 
erties from modern coherence theory are demonstrated to be preserved 
under parametric down-conversion. Based on the generalised Siegert re¬ 
lations and the Cauchy-Schwarz inequality a measure of entanglement for 
the two photon amplitude is proposed. Upper bounds for this measure are 
found in terms of the golden ratio for maximally entangled states. Two 
inequalities are derived for the two photon amplitude for which there exist 
a transition zone from super-Poisson statistics to sub-Poisson statistics for 
down-converted partially coherent fields. 


1 Introduction 

The two photon amplitude (TPA) plays a role of paramount importance in 
the concept of quantum entanglement which has been largely recognized as a 
fundamental property of the non linear phenomena known as parametric down 
conversion (PDC) |T]. The widespread use of bi-photon generation techniques 
and the improvement of coincidence detection of photons has led to the un¬ 
derstanding of the entanglement of bipartite quantum systems as well as to 
the implementation of efficient quantum information protocols. However in the 
vast majority of investigations, the pump field is assumed to be coherent and 
the strategy consists in tailoring the TPA to investigate changes in the spectral 
properties of bipartite systems that can generate high dimensional entanglement 
or lead to new entanglement schemes The aim of the present paper 

is to include the correlation properties of the pump field in order to investigate 
correlation-induced entanglement effects. 

TPA shape manipulation has been investigated by two main approaches, by 
changing the pump spectrum and by adjusting the phase matching conditions. 
The angular spectrum of the pump laser, preserved it the two photon state, 
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allows taking control of the transverse correlations by transferring the image 
from the pump to the down-converted field citePhysRevA.57.3123. Focusing 
a laser pump field has led to the evaluation of the correlation function in an 
imaging system with two photon geometric optics effects by passing the laser 
beam through a focusing lens before the down-conversion crystal ,7j. Down- 
converted photons have been used to investigate the effect of pump spectrum 
and the various parameters of entangled-photon imaging systems on the imaging 
resolution [5]. Entanglement enhancement has been investigated on the basis 
of a judicious selection of eigenmodes in angular momentum entanglement [6|. 
The characterization of entanglement by determining an appropriate set of bi- 
orthogonal mode pairs has been performed by using the Schmidt decomposition 
of the TPA |Jj. This approach has been used for engineering quantum correla¬ 
tions by controlling spectral amplitudes m- The eigenvalues and eigenmodes 
of the density matrix have been used to approximate the TPA by a bi-Gaussian 
function and for calculating the Schmidt number HU- 

The two photon state can be controlled by using transversely patterned 
quasi-phase matching gratings to tailor the spatial mode function that describes 
the state of photon pairs in PDC [12]. It has been shown that type-II down- 
converted photons with broadband pump exhibit a continuous superposition of 
two-photon states in which the probability amplitude for each state depends on 
the pump bandwidth whose envelope function determines the range of pump 
energies and the phase-matching function determines the energy spectrum Ca¬ 
in all the approaches mentioned above, the common characteristic is the 
use of coherent sources from down-converted fields are obtained. An important 
number of potential sources, namely, partially coherent sources have not been 
included and only a few studies have put forward the connection between en¬ 
tanglement and partial coherence. The mathematical theory of duality between 
partial coherence and entanglement comprises the description of propagation, 
diffraction and interference of bi-photon fields in optical linear systems di- 
coherence properties of down-converted fields were investigated in terms of the 
physical properties of a coherent pump beam [15] , In this connection, recently 
it was proved that down-converted photons retain the coherence properties of 
the pump held and the maximum entanglement is bounded by the degree of 
spatial coherence m- With the development of modern coherence theory, in 
the last decade a number of results and techniques have been established that 
brings deep insight into the concepts of coherence and polarisation expressed 
through the cross spectral density matrix [T71 [151 EH HOI [HI HH [HI [HI dl d] . 

In the present paper theoretical implications of correlation-induced entan¬ 
glement effects, i.e., tailoring the TPA photon amplitude with partially coherent 
fields on bipartite partially entangled systems, are investigated. The paper is 
organized as follows. In section 2 it is demonstrated that the coherence proper¬ 
ties of partially coherent pump fields are preserved in the two photon entangled 
fields. In section 3 the cross spectral density operator is introduced. The ex¬ 
pected value of this operator is related to the two photon amplitude trough the 
kernel of the one photon amplitude. In section 4 an extension to the measure of 
entangled states proposed in [26] is established for the TPA of partially coherent 
fields that satisfies the Siegert relation. Under this formalism it is found that 
maximally entangled states have an upper bound given by the golden ratio. An 
important consequence is that entanglement beyond the golden ratio implies the 
violation of the Cauchy-Schwarz inequalities for maximally entangled states and 
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a transition from super-Poisson statistics to sub-Poisson is involved. Finally, in 
section 5 conclusions are presented. 

2 Coupling of partially coherent down-converted 
fields 

The TPA contains all the information of the down converted fields field at the 
moment of birth inside the non-linear material. The correlation time of the 
pump field was believed to play no role in the two photon second order corre¬ 
lations, fundamentally because in Mandel’s original derivation the bandwidth 
was taken to be zero and in part because the concepts of partial coherence in 
the space frequency domain was not fully developed. Note that Wolf’s theory 
of partial coherence in the space frequency domain was published a year later 
than Mandel’s theory, in 1986 mi27!. In the vast majority of investigations 
concerned with two photon systems, the pump beam is assumed to be coherent. 
If the pump is partially coherent the two-photon field operator depends on the 
interaction time, the representation of the classical field and the non linear prop¬ 
erties of the crystal. The expected value of the field operator depends on both 
the expected value of the quantum field operator and the cross spectral den¬ 
sity matrix. Partially coherent fields pertain to a large class of electromagnetic 
sources that obey super-Poisson statistics although not fully unbiased, they are 
often referred to thermal fields. 

Among the large class of partially coherent sources, here we are interested in 
partially coherent fields for which the fourth order correlation function (iq, r 2 ) 
can be expressed through the second order correlation function, i. e. those for 
which the Siegert relation holds [25] 


r( 2 )(r 1 ,r 2 )=rW(r 1 ,r 1 )r«(r 2 ,r 2 ) + |rW(r 1 ,r 1 )| 2 . (1) 

With this relations the fourth order correlation functions are characterised 
in terms of second order correlation functions. 

In the following the down-conversion of photons from the partially coherent 
pump field will be analysed following Mandel’s original derivation pQ within the 
time where the interaction takes place. The analysis is valid for non-degenerate 
type I and type II SPDC. Consider an electromagnetic field Vj(r, t) entering a 
dielectric non linear crystal with bilinear susceptibility y, the Hamiltonian of 
the system is 


H= \ J XijiUi(T,t)UjMV,(r 0 ,t)6r, (2) 

where D is the spatial domain where the interaction of light and the crystal 
takes place and Vj(r 0 ,f) = V;eM ko ' r °~ w ( k °Wl is a member of the classical pump 
field ensemble of realizations with average frequency uiq. The quantized field is 
given by 


U{ T,t) 
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where a 1 (k, t) is the field creation operator of the mode k, ets is a unit polari¬ 
sation vector depending on the wavevector k, s = 1,2 is the polarization index 
and r spans the transverse plane perpendicular to the propagation axis. The 
Hamiltonian is then given by 


H = J D J 

e i[(k 0 -k-k')r-o;(k 0 )t] drdkdk , ( 4 ) 

By defining A = l/L 3 / 2 ae*“* k ^ the Heisenberg equations of motion can be 
written as follows 

A= zKvJ d J j ^oi e ks,k's'Al s ' V i( r o,t) (5) 

e *[(k 0 -k-k')T] e -j[(a;(k 0 )-w(k)-w(k , ))t] drdkdk; / 

where eksi.k's'j represents the unit polarisation vector. In order to determine 
the field mode the Heisenberg equations of motion must be integrated over the 
interaction time At. If the pump field is considered coherent with a coherence 
time of the order of microseconds and since the interaction time is typically 
a fraction of a nanosecond the time integral is straightforward and yields the 
familiar sine function, whose argument contains the phase matching conditions, 
of the frequency and the interaction time. In the case of partially coherent 
sources the situation is different. The values of temporal coherence ranges from 
a few to hundreds of femtoseconds for partially coherent sources and fractions 
of nanoseconds for multimode lasers. With the development of the modern co¬ 
herence theory |29] , an important number of sources with prescribed properties 
of coherence has been put forward so that the range of temporal coherences 
available under the time integration is non negligible. Hence the integration 
over time is expressed as follows 


A= i^LLI / ( 6 ) 

e i[(k 0 —k—k')-r] e —i[(w(k 0 )—w(k)—w(k'))t')] drdkdk ' d ^_ 

It is convenient to express this integral as follows 


A = 


l 

ihLA 



V[(r, dt' 



ij £S e- l[k °- k - k ' lr dr 


XijiTij(ks, kV)dkdk', 


(7) 


where St corresponds the interaction time. If the pump field is partially coherent 
the contributions to the response operator of down-converted light at point r and 
time t resulting from the creation operator will be affected by the correlations 
of the input beam within the non-linear material. To show how this happen, 
consider the operator 
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( 8 ) 


U(r,i) = U free 5 t) + f A ks V(r,u)F(k,s,r,t)dk, 

where U f ree is the field operator in the absence of any interaction, V (r, ui) is 
the Fourier spectrum of the input field within the interaction interval St 



and F( k, s, r, t) is given by 

F(k,s,r,t) y J XijiTijOts^'s^dkdk'. (11) 

It was pointed out in the previous section that, by virtue of the Siegert re¬ 
lation r( 2 )(r!,r 2 ) = r^ 1 ^(r 1 ,r 1 )r( 1 ^(r 2 ,r 2 ) + |r( 1 )(r 1 ,r 1 )| 2 , the fourth order 
correlation function of partially coherent beams is fully described by the second 
order correlation function. From the present analysis is clear that the averaged 
measure of light intensity, at points (ri ,t) and (r 2 ,t) for t > St, is given by 
r( 2 )(ri,r 2 ) = (Ut(ri,t)Ut(r 2 ,i)U(r 2 ,t)U(ri,t)), therefore it is fully charac¬ 
terised by the second order correlation function given by 


r (1) (ri,r 2 ;w) = [ f (^j ts A' s ')(V'*(ri,o;)V - (r 2 ,a;)) 

J k J k' 

F*(k, s, k', s', ri, t)F(k, s, k', s', r 2 ,t) dkdk'. (12) 

The first term in the integrand governs the quantum system of down converted 
photons, it depends on the vacuum expectation value of the boson operators 
during the interaction time and the phase matching conditions. The second 
term, recognised as the cross spectral density function of an electromagnetic 
field, is the correlation function that governs the spatial coherence properties 
of the classical pump field at frequency u>. Therefore the expected value of the 
two photon amplitude contains the information about both the down converted 
photons and the classical correlation properties of the pump field. This result 
reveals the dual world to which down-converted photons belongs, when the 
pump field is partially coherent the contribution to the two photon amplitude 
is much richer than the previous case. We shall see in the next section that 
in order to quantify the effect of field correlations in down-converted photons 
the cross spectral density matrix can be conveniently expressed as a density 
operator allowing the evaluation of the second order correlation function. 

3 The cross spectral density operator 

It is well known that the second order correlation function satisfies the Wolf 
equations and the Helmholtz equation. The Wolf equations have a quantum 
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counterpart whose solution are the TPA; as a consequence the fourth order cor¬ 
relation function also satisfies the Helmholtz equation [25] • This remarkable 
result constitutes a further demonstration of the duality between two photon 
systems and partially coherent fields. These equations posses identical math¬ 
ematical structure in both position and Fourier space, besides both share the 
hallmark of quantum-optical systems namely, linearity. The solutions are there¬ 
fore connected and describe mathematically equivalent physical phenomena. 
The second order coherence function represents the oscillatory behaviour of 
partially coherent light in propagation and diffraction. The two photon am¬ 
plitude represents the probability amplitude of oscillation modes in the system 
of down-converted photons. Therefore the two photon Wolf equations imply 
that the TPA admits a modal representation in terms of the eigenvalues and 
eigenvectors of the second order correlation function. 

Despite of this compatibility further clarification is necessary. The second 
order correlation function and the cross spectral density function are correlation 
functions with analogous mathematical definition, however do not describe the 
same optical phenomena. The Young interference experiment reveals the field 
source correlations at the plane of interference. The mathematical description of 
the interference pattern involves two types of coordinates, i.e., the intensity, in 
terms of coordinates from the plane of interference, is weighted by the degree of 
coherence in terms of coordinates from the source plane. This subtle distinction 
reveals the nature of the coherence functions as source correlation functions. The 
degree of coherence does not provide any information of correlations that may 
occur at the interference plane, where the light is detected. The second order 
correlation function describes correlations of arrivals, i.e., detection correlations 
and the same holds for the fourth order correlation function. In this regard do 
not contain information about correlations at the source plane because they are 
in terms of photon counting statistics. For the sake of comparison it is possible 
to calculate the degree of coherence at the interference plane, the result is that in 
general is different from the source. The same holds in the other case, intensity 
correlations can be measured at the source plane but in general these will not 
be equal to the measurements at the detection plane. 

In the following we proceed to elucidate the connection between the second 
order correlation function and the cross spectral density matrix in the context 
of spontaneous parametric down conversion for partially coherent light. One 
of the most important results in vector coherence theory states that when the 
cross-spectral density (CSD) matrix satisfies square integrability: 



Tr W t (ri,r 2 ;aj)W(ri,r 2 ;w) dridr 2 < oo, 


(13) 


hermitian symmetry: 


W f (ri,r 2 ;w) = W(r 2 ,ri;w), 


(14) 


and is non-negative definite: 



D 


f(xi)W(ri,r 2 ;w)f(x 2 )dridr 2 > 0, 


(15) 
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where f (r) is an arbitrary row matrix, hence the CSD becomes a Hilbert-Schmidt 
operator. This is an important result because the CSD functions can be ex¬ 
panded in Mercer’s series (30| j31j 


W(ri,r 2 ;w) = ^ A„^(ri) 0 $„(r 2 ), (16) 

n 

where coefficients A„ are the eigenvalues and <&„(r) = [<p n ,i{ r)] are the eigen¬ 
vectors of the Fredholm integral equation given by 


$n(ri)W(r 1 ,r 2 ;w)dr 1 dr 2 = A„$„(r 2 ). 


(17) 


ID 


Because of relations (131, (14) and (15) the eigenvalues A ra are real and non¬ 
negative, and the eigenvectors 3?„(r) are orthonormal. Therefore, decomposi¬ 


tion (16) represents the cross spectral density matrix of any partially coherent 


and partially polarised electromagnetic field as a superposition of completely 
coherent completely polarised modes. The vector integral equation repre¬ 
sents a set of two coupled scalar integral equations whose kernel W(ri, r 2 ; w) is 
for the general case unknown. However considerable simplification is obtained 
when the kernel is diagonalised, in such a case the integral equations uncouple 
and can be treated independently. 

Solving the integral equation ,for any particular case is not an easy 
task. For those cases in which this can be achieved, a fundamental property 
of Mercer’s series is that they are absolutely and uniformly convergent. The 
eigenvalue distribution typically can involve a large number of elements, this fact 
makes possible to express the mth-eigenvalue for the ith-component of matrix 


(16 1 in terms of power series. By the central limit theorem, the expansion of 


the mth-eigenvalue is of the form 


m \2i 

WAHD-l)'-^. ( 18 ) 

1=0 

The parameter A,, that satisfies 0 < A* < 1, depends on the ratio between the 
effective width of the source a s and the correlation length a c in a functional 
form determined by the particular type of source. The power series expansion 
represents the contribution, per number of photons, of the mth-mode to the 
coherence of the field. The second order correlation function expressed in terms 
of the density operator p is given by 


(V*{ri;u)Vj(r 2 ;u})) = ^(^(^(ri; w)V,-(r 2 ; w)|</> ? ), (19) 

Q 


where 


P= ^2 p n(a)\a){a\. 

n—0 


( 20 ) 


7 







Therefore the modal expansion determines uniquely the contribution of modes 
per number of photons to the coherence of the field. The eigenvalue expansion in 
the coherent mode representation is related to the coherent state representation 
as follows 


\2k 




ra—0 k—0 


By using the photon number operator, which satisfies a^a\m) = m|m), it follows 
from Eq. (21) that the cross spectral density operator, say jin, is given by 


M a — 


e -[ a ki“kiIn(Ai)] 2 _ 


( 22 ) 


Therefore, the expected value of the cross spectral density function is found to 
be 


Tr{pfiu} = e -K+«^ In (A i ) 2 _ 


(23) 


The density operator p is represented as a superposition of coherent states which 
form an overcomplete basis but they are non orthogonal, the lack of orthogo¬ 
nality of the coherent state representation also holds and is revealed by the fact 
that the basis of the eigenvalues power expansion is also non-orthogonal. From 
Eq. (23) it can be seen that if A —1 the field is completely coherent and consists 


of only one mode for a given number of photons, the corresponding CSD can 
therefore be factorized. If A —> 0 the field consists of a superposition of infinite 
modes weighted by their corresponding eigenvalues and turns out to be incoher¬ 
ent. The larger the number of modes, the smaller number of photons per mode 
so that the corresponding CSD can no longer be factorized. In the more general 
and realistic case, the pump field is partially coherent, the expected value is 
then formed by a superposition of several modes and depending on the strength 
of the eigenvalues, the number of photons per mode contribute with different 
amounts to the field coherence. It is important to note that down-conversion 
efficiency relies on the phase matching conditions, when these are fulfilled and 
plane waves are used as pump, the sine function attains its maximum value. If 


the pump held is partially coherent this situation changes, Eq. (23) shows that 


the expected value of the cross spectral density operator attains its maximum 
value for a larger frequency range centred at the phase matching condition. As 
a consequence maximum efficiency can be attained for partially coherent beams 
for a given mean frequency of the pump beam. 

In figure [l] we compare the contribution of the typical sine function versus the 
expected values of the CSD density operator. It can be seen that for completely 
coherent fields the sine function imposes a natural limit for which the states be¬ 
come rather constrained. If partially coherent light is employed, this restriction 
no longer appears and a number of new states are feasible. In particular for 
incoherent fields, i.e., when A,; —> 0, the two photon amplitude approaches to 
a delta function, the coherence properties can therefore be engineered a priori 
opening new possibilities to attain maximally entangled states. In this regard 






Figure 1: Expected values of the cross spectral density operator for different 
values of A. Cohrent case A,; —» 1 is represented by large-dashed curve. Partially 
coherent case Aj = 0.5 is represented by medium-dashed curve. Incoherent case 
Xi —► 0 is represented by short-dashed curve. For comparison, solid line is the 
typical sine function for the coherent case. In figure n = ko — k — k'. 


a relationship between entanglement, the two photon amplitude and partially 
coherent fields will be investigated in the next section. 

Now some consequences of the inclusion of the correlation properties will be 


shown. On substituting Eq. (231 into (24) the second order correlation function 


is 


r (1) (n, r 2 ; w) = JJ 5kk' r R{pP>ii}F*(k, s, k', s', r x t)F(k, s, k', s', r 2 , f)dkdk'. (24) 

The two photon amplitude now exhibits a more involved phenomenon. For 
t > At it describes the spectral contribution of coherence to the probability of 
detecting one photon at r with frequency ui. 

A CSD is called genuine if it satisfies the non negative definiteness condition 
[32], this is a fundamental feature of all physically realizable correlation func¬ 
tions. However not every correlation function associated with an electromag¬ 
netic field satisfies this condition. A necessary and sufficient condition states 
that a correlation function represents a genuine CSD if and only if it can be 
formally expressed in the form mm 


W(ri,r 2 ;w)= J J p(v i,v 2 )H*(ri, vi; w)H(r 2 , v 2 ;w)dvidv 2 . 


(25) 


This important result in modern coherence theory states that the CSD pertains 
to a class of functions in a reproducing kernel Hilbert space (RKHS) and satisfies 
the reproducing property (17). Since both correlation functions in the integrand 
of Eq. (241 satisfy the non negativity requirement, and because (Al,Aw) = Skk', 
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their product can be regarded as the product of two reproducing kernels in the 
more general Hilbert space H = Hq g) Hs . Functions pertaining to the space 
Hq are functions in the quantum ensemble of states and functions pertaining 
to the space Hs are functions of the statistical ensemble of realizations. By 
the integral representation theorem for reproducing kernels [S3], the product of 
the correlation functions in the integrand of Eq. (241 is a reproducing kernel. 
Therefore, on substituting Eq. (22) into Eq. (24) the function 


r( 1 )(n,r 2 ;w)= //e- ( ^ + ^ )In(Al) V*(k,s,r 1 ,t)F(k,s,r 2 ,t)dk 


(26) 


do satisfy the property (25) for Vi = v 2 = k and indeed becomes a reproduc¬ 


ing kernel. Hence the dual property of the quantum nature of down-conversion 
process and classical correlations are contained in the new reproducing ker¬ 
nel Hilbert space of the two photon amplitude. The two photon amplitude 
T (2) (ri , r 2 ; u>) contains all the information of the classical and quantum system 
through the one photon amplitude T^^ri, r 2 ;w) which describes the joint con¬ 
tribution of the classical correlation functions to the two photon amplitude of 
down converted photons. 


It could be tempting a decomposition of the TPA similar to Eq. (16), how¬ 


ever, since the TPA in general is not hermitian, a Schmidt decomposition [34| 
is possible instead. In the next section it will be shown that such a a decom¬ 
position implies a link between partially coherent fields and non-classical fields, 
i.e fields with super Poisson statistics and fields with sub- Poisson statistics. 


4 Partially correlated fields and quantum entan¬ 
glement 

The results of the previous sections can be exploited to investigate the effect of 
correlation properties of the pump beam on the two photon wave function. In 
coherence theory two relationships between the eigenvalues and the properties 
of the CSD matrix, namely 


/'Tr[W t (r 1 ,r 2 ;a;)W t (r 1 ,r 2 ;a;)]dr 1 dr 2 =^A 2 , (27) 

** n 

and 


f f Tr[W t (r,r;w)]dr = ^ A„, 

n 

lead to the definition of the effective degree of coherence 135, 36] 


(28) 


^eff 


En A n 

(E„A„) 2 ' 


(29) 


This value represents a relative measure of coherence within a specific area. 
Note that the term 1/ E) A n corresponds to the Schmidt number which is a 
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measure of entanglement. The effective degree of coherence and the Schmidt 
number seem to play equivalent role for coherence and entanglement. In the 


present case, the effective degree of coherence (29) can be calculated from the 


eigenvalue power expansion (18), substituting Eq. (18) into (29) the effective 


degree of coherence is given by 


faff — 


Er=o(E;= 0 (-i) 


\ 2m \ 2 

-11 m - - ) 

2m!/ 


E“„ELo(-i) 


A 2m \ 2 ’ 

-1 - 


(30) 


2m! / 


In general the eigenvalues depend on the ratio A = and are determined 

fundamentally by the analytic form of the CSD matrix. It is clear that for a 
given genuine CSD this result shows the fact that there exist an effective number 


of modes with a defined number of photons per mode as stated by Eq. (23). 


It may be tempting to represent the TPA in the form (161, however the lack of 


hermiticity imposes a rather strong condition for the implementation of a modal 
expansion therefore the two photon system does not admit a coherent mode 
representation but a Schmidt decomposition. The main distinction between 
these decompositions relies in the fact that the two photon system Hilbert space 
comprises the tensor product of Hilbert spaces of functions u, v € Hq and 
x,y £ Hs- A pure bipartite system is given by 


£ = ax ig) u + (3y <g> v. 

In the Schmidt basis the density matrix of the system is given by |34j 


(31) 


P = ^ = ^2 CiC*aUiU 


) ViV. 


t 


(32) 


which can be expressed as the product of reduced density matrices of the form 


(161. This implies that whereas non hermiticity of operators and non orthogo¬ 
nality of bipartite system basis do not allow us to express Eq. 0 in the form 
(16), the reduced density matrix admits a Mercer’s series provided hermiticity, 


orthogonality and square integrability is fulfilled in each subsystem. An expan¬ 


sion of the form (32) for entangled states of bipartite systems has been proposed 


to express the system as the sum of a completely entangled states and a factor- 
izable state [26| . In the same spirit and by taking into account that for partially 
coherent fields the Siegert relations holds, we propose a Schmidt decomposition 
for the TPA as the sum of a factorizable non entangled amplitude function plus 
a completely entangled amplitude function of the form 


4 2) (r 1; r 2 ) = V^T^n.ra) + yjl - M% T}(n, r 2 ), (33) 

with the entangled component 


Tg(ri, r 2 ) = |r 1 (ri,r 2 )| 2 , 


(34) 
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and the factorized component 


r}(ri,r 2 ) = r 1 (r 1 ,r 1 )r 1 (r 2 ,r 2 ). (35) 

A measure of entanglement that characterizes such a bipartite system is then 
given by M E £ C that satisfies 0 < \M E \ < 1. The corresponding second order 
correlation functions are related by the Cauchy-Schwarz inequality 

-\/Me|r (1) (r,,r 2 )| 2 < v'l-M e 2 rW(r 1 ,r 1 )r«(r 2 ,r 2 ). (36) 

In order that this inequality holds, the inequality y/M e < y/l — must also 
hold, therefore the degree of entanglement is bounded between the interval 


0 < M e < 


(-1)±V5 


(37) 


There are two values for which the inequality holds, one real and one complex, 

1 -j- y/5 

interestingly they are in terms of the golden ratio <f> = —-—. Since <j> — 1 = 


y/E — 1 . 


in the real case the degree of entanglement satisfies 


0 < M e < (38) 

9 

the negative solution, associated with the complex plane, implies that the com¬ 
plex degree of entanglement, denoted by M e , satisfies 


0 < M e < i<j>. 


(39) 


Beyond this bound, the Cauchy-Schwarz inequality no longer holds, conse¬ 
quently the TPA may undergo changes in its statistics due to entanglement 
effects. This remarkable result indicates that the degree of entanglement, within 
its interval of definition, can lead to changes in bi-photon counting statistics. 
It is interesting not only because there exists an upper bound for entanglement 
but because establishes a threshold between entanglement and the violation of 
the Cauchy-Schwarz inequality which is a key point for the anti-bunching effect. 
The nature of entanglement can be enlightened further by characterising the 
TPA in terms of the second order correlation functions. On substituting Eqs. 


(33), (34) and (35) into Eq. (36) the following pair of inequalities are found 


2 V / Me|r (1) (ri,r 2 )| 2 < r( 2 )(n,r 2 ) < 27T^Mfr«(r 1 ,r 2 ). (40) 

The violation of the Cauchy-Schwarz inequality is a necessary but not sufficient 
condition for anti-bunching effects hence M e = f does not guarantee anti¬ 
bunching. Specifically the sufficient condition implies that T01(r 1; r 2 ) < 1, 
beyond this bound the statistics turns out to be sub-Poisson with a pump field 
originally super-Poisson or Poisson in the coherent case. The value for which 
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Degree of entanglement 


Figure 2: Proportion between the maximal entanglement (blue line) and mini¬ 
mal entanglement (purple line). Region limited by parallel red lines represents 
the transition zone, where Cauchy-Schwarz inequalities are violated, from super- 
Poisson statistics for M e < <j> -1 to sub-Poisson statistics for M e > . 


this occurs is found to be M e > V?>/2. Figure [ 2 ] illustrates the golden section 
<j> < M e < \/3/2 where this transition phenomena arise. 

This result states that, provided entanglement and coherence properties are 
properly engineered, anti-bunching effects with partially coherent light are feasi¬ 
ble since the fourth order correlation function admits a Schmidt decomposition 
that can be carried out by means of the coherent-mode representation of the 
second order correlation functions. It transpires that coherence properties of 
electromagnetic fields can be used for tailoring the two photon amplitude for 
quantum entanglement strategies and vice versa, quantum entanglement can be 
used to manipulate the photon counting statistics of partially coherent down- 
converted fields. In addition, the obtained result suggests that the non-genuine 
correlation functions, those for which the non negativity condition (15) is not 
satisfied, can be studied in connection with non-classical fields, i.e., those for 
which anti-bunching effects occur, provided the appropriate statistics for pho¬ 
ton counting events is selected. The presented results indicate that correlation 
properties of the pump field introduces a new degree of freedom for tailoring 
the two photon amplitude in entanglement applications, end moreover entan¬ 
glement can be engineered to test transition states between super-Poisson to 
sub-Poisson statistics. 


5 Conclusions 

We have shown that the coherence properties of the electromagnetic pump field 
are preserved under parametric down-conversion. It was demonstrated that 
the two photon amplitude contains full information about both quantum down- 
converted photons and the classical correlation properties of the pump beam. 
By using the coherent mode representation, the statistical nature of the CSD 
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matrix has been transferred to a CSD operator whose expectation value allows 
us to determine the shape of the TPA for a partially coherent pump field. The 
relationship between the expected value and the kernel properties of the CSD 
led to the deduction of the reproducing kernel property of the second order 
correlation function. This property and the modal decomposition account for a 
complete characterization of the TPA in terms of coherent modes. As a conse¬ 
quence, by choosing appropriately an orthonormal basis of the pump beam, the 
TPA can be shaped in terms of the number of modes resulting in a new strategy 
for entanglement engineering. We put forward a measure of entanglement by 
means of the TPA as a direct extension of the entanglement measure based on 
states This measure allows us to establish upper bounds for the degree of 
entanglement, the most important conclusion is that field correlations through 
the modal decomposition can be deployed to tailoring the TPA and conversely 
the degree of entanglement can be used to test the transition of fields obey¬ 
ing super-Poisson statistics to sub-Poisson, therefore anti-bunching effects are 
expected provided the appropriate correlation properties are engineered. The 
presented results indicates that correlation properties of the pump field intro¬ 
duce a new degree of freedom for down-converted fields. This results can be 
useful in quantum communications and quantum information theory. 
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